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Abstract 

The evolution equation for the shear is reobtained for a spherically 
symmetric anisotropic, viscous dissipative fluid distribution, which al- 
lows us to investigate conditions for the stability of the shear-free 
condition. The specific case of geodesic fluids is considered in detail, 
showing that the shear-free condition, in this particular case, may be 
unstable, the departure from the shear-free condition being controlled 
by the expansion scalar and a single scalar function defined in terms 
of the anisotropy of the pressure, the shear viscosity and the Weyl 
tensor or, alternatively, in terms of the anisotropy of the pressure, the 
dissipative variables and the energy density inhomogeneity. 
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1 Introduction 



The relevance of the shear tensor in the evolution of selfgravitating systems 
and the consequences emerging from its vanishing has been brought out by 
many authors (see [l]-[6] and references therein). Particular attention de- 
serves the possible role of shear in the violation of the cosmic censorship, 
leading to the appearance of a naked singularity. This has the important 
implication that the resultant spacetime singularity of collapse could become 
visible to faraway observers in the universe (see [4] for a detailed discus- 
sion). Furthermore as it has been recently shown [7] the shear-free flow (in 
the nondissipative case) appears to be equivalent to the well known homol- 
ogous evolution. It should be recalled that homology conditions are of great 
relevance in astrophysics [8]-[10]. 

Accordingly, it is quite pertinent to ask under which conditions the shear- 
free condition is stable. This question in fact entails two different (but re- 
lated) questions, namely: 

• Under which conditions an initially shear-free flow remains shear-free 
all along the evolution ? 

• Under which conditions a fluid with a small initial shear evolves, keep- 
ing always a small value of shear (Lyapunov stability)? 

Now, in the study of self-gravitating compact objects it is usually as- 
sumed that deviations from spherical symmetry are likely to be incidental 
rather than basic features of the process involved. Thus, since the semi- 
nal paper by Oppenheimer and Snyder [11], most of the work dedicated to 
the problem of general relativistic gravitational collapse, deal with spheri- 
cally symmetric fluid distribution. Accordingly we shall consider spherically 
symmetric fluid distributions. 

The purpose of this work is to provide answers (at least partial) to the 
above questions. Obviously for doing so we need an evolution equation for 
the shear. Such an equation was derived by Ellis [12], [13], for a perfect fluid 
without any kind of symmetry. Here we shall reobtain that equation for 
a spherically symmetric anisotropic viscous dissipative fluid. The physical 
motivation to consider such a general fluid has been explained in detail in 
[14]-[17] (and references therein). However, since the study of the general 
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case is quite complicated and will probably require the use of numerical 
techniques, we shall consider in detail the case of a geodesic fluid. 

Nevertheless, even though the restriction to the geodesic case is dictated 
by the sake of simplicity in the analysis, the physical relevance of such a case 
becomes evident when we recall that Friedman-Lemaitre-Robertson- Walker 
(FLRW) metrics are shear-free. In other words, it appears that for a geodesic 
fluid (without dissipative fluxes), the stability of the shear-free condition is 
somehow equivalent to the stability of FLRW spacetime itself. 

For the geodesic case it will be shown that departures from the shear-free 
condition are controlled by the expansion scalar and a scalar function, which 
is defined in terms of the the Weyl tensor, the anisotropy of pressure and the 
shear viscosity or, alternatively, in terms of pure physical variables. Such a 
scalar function appears in a natural way in the orthogonal splitting of the 
Riemann tensor and is related to the Tolman mass [17]. Together with the 
expansion scalar, this scalar function controls the evolution of the shear. 

In the next section we shall present the general equations and definitions 
to obtain the evolution equation for the shear. Next in Sec. Ill that evolution 
equation is obtained. The origin and some properties of the scalar function 
which plays such an important role in the evolution of the shear are discussed 
in Sec. IV, and in Sec. V the obtained equation is used to analyze the geodesic 
case. Finally, a discussion of results is presented in the last section. 

2 ENERGY-MOMENTUM TENSOR, KINE- 
MATICAL VARIABLES AND FIELD EQUA- 
TIONS 

We consider a spherically symmetric distribution of collapsing fluid, bounded 
by a spherical surface S. The fluid is assumed to be locally anisotropic (prin- 
cipal stresses unequal) and undergoing dissipation in the form of heat flow (to 
model dissipation in the diffusion approximation), null radiation (to model 
dissipation in the free streaming approximation) and shearing viscosity. 

Choosing comoving coordinates inside S, the general interior metric can 
be written 

ds 2 = -A 2 dt 2 + B 2 dr 2 + R 2 (d6 2 + sin 2 ## 2 ), (1) 
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where A, B and R are functions of t and r and are assumed positive. We 
number the coordinates x° = t, x l = r, x 2 = 9 and x 3 = <p. 
The matter energy-momentum T a p inside X has the form 

T a/ 3 = (/J,+P±)V a Vp+Pj_g a p+(P r -P ± )xaXp+qaVp+V a qp+elJp-2r}<j a p, (2) 

where /i is the energy density, P r the radial pressure, Pj_ the tangential 
pressure, q a the heat flux, e the energy density of the null fluid describing 
dissipation in the free streaming approximation, r\ the shear viscosity coeffi- 
cient, V a the four velocity of the fluid, x a a un h four vector along the radial 
direction and l a a radial null four vector. These quantities satisfy 

V a V a = -l, V a q a = 0, X a Xa = l, X a V a = 0, l a V a = -l, l a l a — 0. (3) 

Observe that we have assumed the shear viscosity tensor 7r a/3 to satisfy 
the relation 

7r a/3 = -2r]a af} , (4) 

where a a/ 3 is the shear tensor. However this last equation is valid only within 
the context of the standard irreversible thermodynamics (see [18], [19] for 
details). 

In a full causal picture of dissipative variables we should not assume (4). 
Instead, we should use the transport equation derived from the corresponding 
theory (e.g. the Muller-Israel-Stewart theory [20]-[22]). However for the 
sake of simplicity, in this work we shall restrict ourselves to the standard 
irreversible thermodynamics theory. 

The acceleration a a and the expansion O of the fluid are given by 

a a = V a . fi V f> , Q = V a ;a . (5) 

and its shear cr a/ g by 

o"a/3 = V( Q;/9 ) + a( a Vp) - -Qh a p, (6) 

where h a p = g a p + V a Vp. 

We do not explicitly add bulk viscosity to the system because it can be 
absorbed into the radial and tangential pressures, P r and P±, of the collapsing 
fluid [23]. 
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Since we assumed the metric (1) comoving then 

V a = A" 1 ^, q a = qB^Sf, l a = A" 1 ^ + B- l 5^ X a = B^S?, (7) 

where q is a function of t and r satisfying q a = q\ a - 

From (5) with (7) we have for the acceleration and its scalar a, 

A' ( A' \ 2 

ai = -J, « 2 = a a a a = [j^J > (8) 

where a a = a% a , and for the expansion 

where the prime stands for r differentiation and the dot stands for differen- 
tiation with respect to t. With (7) we obtain for the shear (6) its non zero 
components 

g 11 = \b 2 <t, a 22 = ^§- = -\R 2 a, (10) 
3 sm 6 3 

and its scalar 

a"*** = \°\ (11) 

where 

1 (B R\ 

°=a{b-r)- (12 > 

Then, the shear tensor can be written as 

<7aP = ° (xaX/3 ~ \ h <*p) ■ ( 13 ) 

2.1 The Einstein equations 

Einstein's field equations for the metric (1) are given by 

G af} = 8nT a(3 , (14) 
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its non zero components with (1), (2) and (7) become 

8nT 00 = 8n((i + e)A 2 
R" 



2— —I - 

B + RIR 



8nT< 



01 



A\ 2 

BJ ~ R 

8n{q + e)AB 



R 



B R \RJ 



8nT, 



R! BR! RA'\ 
B~R ~ ~RA) ' 
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8nT 22 = -^-^33 = 8tt (p ± + - V a) R 2 
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B R A (B K 
B + R~ A\B + R 
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~A + ~R 
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I fli? 
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(15) 
(16) 



(17) 



(18) 



2.2 The mass function 

Let us now introduce the mass function m(t,r) [24] (see also [25]), defined 
by 



p3 

m — -^-^23 — — 



R 
~2 



R/ 



(19) 



Following Misner and Sharp [24], It is useful to define the proper time 
derivative Dt given by 

1 d 

Dt = ~, (20) 



and the proper radial derivative D R , 



D 



R 



Adt 



l_d_ 



(21) 
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where R defines the areal radius of a spherical surface inside S ( as measured 
from its area). 

Using (20) we can define the velocity U of the collapsing fluid as the 
variation of the areal radius with respect to proper time, i.e. 

U = D T R < (in the case of collapse). (22) 

Then (19) can be rewritten as 

„ R' /„ tt2 2m\ 1/2 

E= — = (l + U 2 -—) . (23) 



B V R ) 

Using (15)-(17) with (20) and (21) we obtain from (19) 



D T m = —An 

and 



P r -i w )U + ( 1 E 



R 2 , (24) 



D R m = 4n(ji + q^R 2 , (25) 



which implies 



rK / (J\ 

m = Air J [Ji + q—j R 2 dR, (26) 

(assuming a regular centre to the distribution, so m(0) = 0). Integrating 
(26) we find 

3 THE EVOLUTION EQUATION FOR THE 
SHEAR 

We shall now proceed to deduce the Ellis evolution equation for the shear, 
for the specific fluid distribution discussed in the previous section. For that 
purpose it will be convenient to express the energy momentum tensor (2) in 
the equivalent form 



-a/3 



pV a Vf, + Ph aP + n Q/3 + q (V aX p + XaVp) - 2 Wa p (28) 
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with 

p= Pr + 2P± 



3 

ft — + e, 

9 = Q + e, 
n = P r - P ± , 

n a/3 = n (xaXp - l^Kp 

Now, Ricci identities for the vector V a read 

K,/3v V H = V <*$;v ~ V a]V ./3, (29) 

or, using the well known expression (remember that vorticity vanishes due 
to the spherical symmetry) 

V a ;n = -a a V^ + o aii + -6/i aA4 , (30) 

we obtain 

2 R aP» V p = a <*;[P V tA + a °y^-A + ^[/Sim] + + ^©^[/M- (31) 

Contracting Eq.(31) with V^g afl , we find the Raychaudhuri equation for 
the evolution of the expansion 



®; a V a + i© 2 + a^a af} - a? a = -Vy% = -4vr(A + 3P). (32) 



On the other hand, contracting (31) with V^h^h^ we have 

- KVfJa^ + \eh lf} ), (33) 
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which by using (30), can be written as 



a 



"71/ 



y 



(34) 



Next, the Riemann tensor may be expressed through the Weyl tensor 
C^pp, the Ricci tensor R a p and the scalar curvature 1Z, as: 

-\R p ^-\n{5 p pg a „-g a p5$. (35) 
Contracting (35) withV^V^/i^/i^ and using Einstein equations, we find: 

VpV^h^KR^ = £ 7 „ + y (/} + 3P) ^ - 47rn 7 , + 8^^,, (36) 

where E 1U denotes the "electric" part of the Weyl tensor (in the spherically 
symmetric case the "magnetic" part of the Weyl tensor vanishes), and is 
defined by 

E a p = C^VV, (37) 
whose non trivial components are 



E 
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E22 
E33 



where 
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A R, 
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(39) 
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Observe that the electric part of Weyl tensor, may be written as: 

E a p = S(xaXp-^h afi ). (40) 
Using (32), (36) can be written as 

(41) 

then, from Eq.(34) and Eq.(41) it follows 
E 1V — AtiH 1u + 87rrja yu = 

h°K (<W ~ V p a a „, p ) - ^/i 7 „ + a,,a u - X -a lv (20 + a) (42) 
or, using (13) 

E 1V — AnIi lu + 8nrja.y U = 

h°h$ (a a;M - V p a aK pj - + a 7 a„ - -<v6 + -cr 2 h lu - c^a 7( g. 

(43) 



Finally, contracting (43) with x l X v we obtain 



S - 4vrn + 87r?7(7 = a f - a* + a 2 - °— - JOa - a-^- , (44) 

3 3 RB 

with /t = / >aX ° and /* = /, a y«. 
Using (7), (44) takes the form 



2 



2„ # 



a' a a 2 2^ 



F TF = 5 - 47TI1 + 8nr,a = ^ - - + a 2 - - - ^Qa - . (45) 

This is the evolution equation for the shear we were looking for (we recall that 
there is only one independent component of the shear tensor in our case). It 
is equivalent (in comoving coordinates) to Eq.(71) in [15] or Eq.(lOl) in [17] 
(notice that E in those references equals — £ here). 

We shall use the above equation to study the stability of the shear-free 
condition for the case of a geodesic fluid. However before doing that we shall 
discuss about the origin and the physical properties of Ytf- 
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4 ON THE ORIGIN AND PROPERTIES OF 

Y TF 

As we mention in the Introduction, the scalar function Y TF appears in a 
natural way in the orthogonal splitting of the Riemann tensor (see [17] for 
details). 

Indeed, the electric part of the Riemann tensor (which is one of the ele- 
ment of that splitting) is defined by [26], [27] 

Yaf, = RonpsVW 6 , (46) 

which in turn can be splitted in terms of its trace and the corresponding 
trace-free tensor, i.e. 

Y af , = ±TrYh af) + Y <a f i> , (47) 

with TrY = Y« and, 

/ TrY \ 
Y <ttfi> = Kffp (Y, u - — Vj , (48) 

this last tensor may also be written as 

Y<af3> = YTF (xaXp ~ \ h *p) ■ ( 49 ) 

To obtain TrY and Ytf we may proceed as in [17] or, directly from (46)- 
(49), using the expressions for the Riemann tensor components in terms of 
the Einstein tensor components, given in the Appendix of [16]. Either way 
the result is 

TrY = Y T = 4tt (p, + 3P r - 2Il) , (50) 

and 

Y TF = 4:7r(-U + 2ria)+S. (51) 
Next, using (15), (17), (18) with (19) and (39) we obtain 

»7T7 

— = An^-U + 2 V a)-S, (52) 
which combined with (27) and (51) produces 
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Y TF = -SttII + lQn W + ^ j* R 3 (d r ~h - 3q^j dR. (53) 

Thus the scalar Ytf m &y be expressed through the Weyl tensor and the 
anisotropy of pressure or in terms of the anisotropy of pressure, the density 
inhomogeneity and the dissipative variables. It is worth recalling that a link 
between Y TF and the Tolman mass has been established in [17]. 

We shall now bring out the role of this scalar function in the stability of 
the shear-free condition for the geodesic fluid. 



5 THE GEODESIC FLUID 

If we assume the fluid to be geodesic (i.e. a M = = a) then it follows from 
(8) and rescaling t, that A — 1. In this case (45) reads 

* + Y + ^t + Ytf = ' (54) 

where Ytf is defined through (51) or (53). 

Let us first assume that Ytf — 0, then the general solution of (54) is 

a(r ' t] = [c(r)t - ell edt (l - §c(r) /e*e-§/ edt ~d*)] ' ^ 

where c(r) is a function of integration. After some simple manipulations, 
this solution can be rewritten as 

a(r,t) = / \ r — -~ • 56 

c( r )fe-Wf edt dt-l 

If we now demand the shear to vanish initially, (i.e. cr(0, r) = 0) then as 
it follows from (55) or (56) we must have (assuming 6 is a regular function 
of its arguments) 

c(r) = 0, (57) 

implying a = for all t. Therefore, if Ytf = 0, the only solution compatible 
with an initially shear-free flow is a shear-free flow. This result is also evident 
from (54). 



12 



Before proceeding further, the following remark is in order: Observe that 
the condition Ytf = does not imply conformal flatness (£ = 0) as it is 
obvious from (51), unless we assume further the fluid to be perfect. In this 
latter case the result above is compatible with the one obtained by Stephani 
[28], [29] stating that the most general conformally flat perfect fluid solution 
differs from FLRW only by having nonzero acceleration. 

Let us now consider the case (always with Y TF = 0) when, initially, the 
flow is close to the shear-free condition, but still a ^ 0. Then the question is 
under which conditions the fluid will evolve, keeping close to the shear-free 
regime (Lyapunov stability). 

Thus, let us assume 

a(0,r) = ea(r), |e| « 1 (58) 

Now, for sufficiently small e and assuming to be a regular function of 
t, we obtain from (55) 

c(r) « 0(e), (59) 

implying 

a^-ee~y &dt . (60) 

From the above it is evident that if > then the shear will remain 
always close to the initial ("quasi-shear-free" condition). However it is also 
evident that for < the situation radically changes and in principle de- 
partures from the quasi-shear-free condition may be expected, depending on 
0. 

Let us now consider the case Y T f ^ 0. 

We shall first assume that Y T f (at least initially) is small, i.e. 

Y TF = aY TF ; a « 1, (61) 

and 

a(t,r) = a (t,r)+pa 1 (t,r), (3 « 1 (62) 

e(t,r) = e o (t,r)+te 1 (t,r), £«1 (63) 

where a (t,r) and o (t,r) correspond to the general solution of (54) for 
Ytf = given by (55). 
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Introducing (62) and (63) into (54) and linearizing we obtain for o~i(t,r) 
the following equation 

o-i + \ (<?o + e ) a x + Y^ ia ° + ]j ?TF = °' (64) 

which may be easily integrated to obtain 
ai (t,r) = /( r ) e -( 2 /3)J> +e )<i t 

_ e -(2/ 3 ) / (a 0+ e o)d t J 1 (j£ 01 a o + a? TF ) e ( 2 / 3 ) / ( CT °+ 6 ^, 

(65) 

where /(r) is an integration function. Introducing this last equation into 
(62) we have 

a(t,r) = a (t,r) + pf( r ) e -W)S(°o+e )dt 

_ e - (2 /3) / (* 0+ e o)d t J (J^e^ + e^ 3 ) / ^ +e ^ di dt. 

(66) 

Let us now assume that the fluid is initially shear-free, which implies 
cr = 0. Then it follows at once from (66) 

a(t,r) = Pf(r)e-^f @ ° dt 

_ e -(2/3)/e dt J a Y TF eW$ e ° di dt. 

(67) 

If, Ytf = (a = 0) then from the condition of an initially shear-free fluid, 
we have that f(r) = 0, implying a(t, r) = for all t, which is the result 
obtained before. 

However, for any a ^ 0, it appears from (66) that the fluid may evolve to 
a non-vanishing shear regime, for any sign of 0, depending on the specific 
time dependence of Ytf- hi other words, even for small values of a the fluid 
may deviate from its initial shear-free condition. 

If the fluid initially satisfies the quasi-shear-free condition, which implies 
that ctq is small but non-vanishing. Then in the case Ytf = a = we obtain 
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as before that for G < departures from the quasi-shear-free condition may 
be expected, depending on 0. 

In the general (a ^ 0) case however, departures from the quasi-shear-free 
condition may be expected along the evolution, depending on I^f for any 
sign of 0. 

An alternative (and useful) expression for the shear, which allows to study 
the same problem, may be obtained as follows. 
From the geodesic condition 

a = A = 1, (68) 
the field equation (16) becomes 

R' B R' _ _ . „ „ . 



or, 

B R! FR 



B R' R' ' 

where F = AnqB. 

Then, (70) with (9) and (12) produces 



Introducing (71) into (54) we get 

a+[(ln(R'R))-™ 
which after integration yields 

J (FR/R')dt 



eJ 

a = 



(70) 



2e+ ^ 3 (l + i-^)^ 3(MM)) - 3 ^ (71) 



a + Y TF = 0, (72) 



J R'RY TF e- 1 ( FR / R ') di dt + C(r)) . (73) 



RR> 

In the non-dissipative case (F — 0), (73) takes the very simple form 



1 

RR 



- t {- j R'RY TF dt + C(r)y (74) 
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Now, if we assume that Ytf = then it follows at once from (73) that 
for an initially shear-free fluid we must have C(r) =0 implying a(t, r) = 0. 

Also, for an initially quasi-shear-free fluid (with Y TF = 0) it should be 
clear that departures from that initial condition are possible if < and 
such that it produces a sharp decreasing of R (remember that in order to 
avoid shell crossing singularities we must have R' > 0). Thus we recover our 
previous result. 

In the general case (Y T f ^ 0) it is obvious from (73) that departures from 
either the shear-free or the quasi-shear-free initial conditions are possible. 



6 DISCUSSION 

We have carried out a study on the stability of the shear-free condition based 
on the evolution equation for the shear presented in Sec. Ill . As shown in the 
previous section, a major role in such study is played by the scalar Ytf- This 
scalar, which appears in the orthogonal splitting of the Riemann tensor, is 
not only related with the Tolman mass as shown in [17], but may be expressed 
through purely physical quantities as in (53). 

We shall not insist on the relevance of the question considered in this 
work (the stability of the shear-free condition) since that was clearly stated 
in the Introduction. However we would like to comment further on the results 
concerning the geodesic case. 

Indeed, as mentioned before, it is worth noticing that in the non-dissipative, 
isotropic case, it may be easily shown that (e.g. see [15]) 

Y TF = £ = a = 0. (75) 

Thus for that particular case (geodesic, non-dissipative, isotropic pres- 
sure) the shear makes the difference between FLRW metrics (Y TF = a = £ = 
0) and LTB metrics (Y TF ^ 0, a ^ 0). 

Therefore, the discussion in Sec.V illustrates for this particular case, how 
departures from an initial FLRW spacetime are controlled by Ytf- 

Finally, it should be emphasized the fact that, even though it may be 
intuitively obvious that dissipative processes, local anisotropy of pressure 
and energy density inhomogeneity should affect the stability of the shear- 
free condition, it is not so obvious that the above mentioned factors affect the 
stability of the shear-free condition only through their specific combination 
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given by (53). Thus for example we could consider the peculiar case when 
all those factors are present but they cancel each other in (53), producing 
Y TF = 0, in such a case the shear-free condition would be stable in spite of the 
fact that the fluid is non-homogeneous, non-isotropic and non-dissipative. 
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